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We derive constraints on the statistics of the charge transfer between two conductors in the model 
of arbitrary time-dependent instant scattering of non-interacting fermions at zero temperature. The 
constraints are formulated in terms of analytic properties of the generating function: its zeroes 
must lie on the negative real axis. This result generalizes existing studies for scattering by a time- 
independent scatterer under time-dependent bias voltage. 



1. Introduction. The concept of full counting statis- 
tics (FCS) for charge transfer in coherent conductors has 
been introduced in Ref. While the average current 
through small junctions may usually be understood in 
classical terms, its fluctuations reflect the quantum na- 
ture of electrons; even in the model of non-interacting 
electrons, the fluctuations of the current are non-classical 
due to the electronic Fermi statistics [2]. Generally, one 
studies the statistical distribution of the charge transfer 
between several conductors coupled by a time-dependent 
scattering matrix. In the simplest case (considered in 
the present paper), there are only two leads, and there- 
fore the charge transfer may be statistically described by 
the probabilities pq to transfer exactly q electrons from 
the left to the right lead (here q is an integer number, 
either positive or negative). The leads are assumed to 
provide a thermal source of incoming electrons with the 
Fermi occupation number np[E). In the context of co- 
herent quantum manipulation, one of the most interest- 
ing setups is that of an "adiabatic pumping" where the 
scattering matrix varies slowly compared to the charac- 
teristic scattering time 0\ . Neglecting the scattering time 
amounts to considering the scattering matrix to be inde- 
pendent of the energy E, i.e., local in time |4i]. In this ap- 
proximation, in the model of non-interacting electrons, a 
general expression (the so-called "determinant formula" ) 
for the charge-transfer statistics {pq} in terms of the time 
evolution of the single-particle scattering matrix S{t) has 
been derived P, [HI . 

Remarkably, a recent progress has been made in under- 
standing the implications of the old results for the FCS in 
the particular case of the bias voltage applied to a fixed 
scatterer Q . In this "bias- voltage" problem, the scatter- 
ing matrix S{t) is changed along a one-dimensional tra- 
jectory parameterized by the U(l) phase (f){t) = / V{t)dt, 
where V{t) is the applied voltage. Within this model, 
non-trivial constraints have been derived for the charge- 
transfer statistics. Namely, not every set of probabilities 
{pq] is allowed, but only those for which the zeroes of 
the generating function 

X(A) = Pi^'^'^ (1) 

g— — oo 



obey certain restrictions (see our discussion below). 

2. Main result. Inspired by this work |6|, we derive 
similar (weaker) constraints for the more general case 
of adiabatic pumping [3]: an arbitrary time-dependent 
scattering matrix S{t) with any number of channels. We 
restrict our consideration to the case of two conducting 
leads and to zero temperature. Under those conditions, 
we find that the generating function ([T]) may take zero 
values only for a discrete set of A such that u = e^^ 
belongs to the negative real axis (— cx), 0). Furthermore, 
we find that the generating function x(A) (and therefore 
the full charge-transfer statistics) is uniquely determined 
by its zeroes Uk, up to an overall integer charge transfer 
x(w) 1-^ M^x(w). These statements constitute the main 
finding of our paper. 

As a simple illustration of the derived constraint, con- 
sider a charge-transfer statistics with the only non-zero 
probabilities p_i^o.i = 1/3- While this statistics appears 
a priori reasonable, it follows from our results that it 
can never be produced for any time-dependent scatter- 
ing matrix S{t) with any number of channels between 
two conductors (since such a statistics would correspond 
to the pair of complex roots u = e'^^'^'^l'^). 

3. Determinant formula and its regularization. The 
generating function ([1]), in the approximation of non- 
interacting electrons and of instant scattering, is given 
by the determinant formula [H, [3, [E [3] ' 

X(A) = det (1 + npiS^e'^P'Se-'^^' - 1)) , (2) 

where S is the time-dependent single-particle scattering 
matrix. At each moment of time, S{t) G U(2M), where 
M is the number of conducting channels in each of the 
two conductors. The operator P3 — diag(lM,OM) is the 
projector counting the charge in the right conductor. The 
operator n^? is the Fermi distribution function in the con- 
ductors. In this paper, we consider only the case of zero 
temperature, and so npiE) = 9{—E) is the step function 
in the frequency representation (in the time representa- 
tion, np{t,t') — i/[27T{t ~ t' + iO)]). Since the electrons 
are non-interacting, we treat them as spinless fermions 
(the spin may be trivially included as a channel index). 

The determinant ([2]) is understood as that of an 
infinite-dimensional operator acting both in the time do- 
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main and in the 2M-channel space. In its original form 
([2]), the determinant is iU-defined: while the operator 
tends to 1 at infinite positive energies, it goes to another 
unitary matrix, S^e'^^^^ Se~^^^^ , at infinite negative en- 
ergies. Therefore the determinant requires a proper reg- 
ularization 0] . A consistent way to regularize the deter- 
minant has been proposed in Ref. 0] (their regularization 
is also equivalent to that used in Ref. 0). At zero tem- 
perature, the Fermi distribution operator is a projector 
n% = np, and the regularized determinant may be writ- 
ten in a particularly simple form 8] , 



X(A) = det {e 



[l-P3 + ^'3^e'^"^^^]}. (3) 



There is a subtle point in calculating the FCS in the 
situation where charge counting is performed over a finite 
time interval: in this case, there are contributions to the 
noise arising from the starting and ending of the mea- 
surement, even in the absence of scattering 

This 

noise grows logarithmically with the observation time. 
To eliminate this contribution from the switching-on and 
switching-off, one often considers a periodic signal S(t) 
(e.g., by repeating the same signal with a large time pe- 
riod). In this periodic formulation, at large number of 
eriods Np, the asymptotic behavior of x(A) is given by 



X(A) ^ [xoiX)] 



(4) 



where xoW = exp[limArp^oo i^lnx(A)] is the counting 
statistics per period. The generating function x(A) may 
be computed with the same formula ([3]), but with the 
finite time interval closed into a loop with the periodic 
boundary conditions. We further consider this periodic 
setup and omit the subscript by using the notation x(A) 
for the charge-transfer statistics per period. 

To write the operators in ([3]) in the matrix form, one 
can perform a Fourier transform in time (introducing dis- 
crete quasi-energies) . In the quasi-energy basis, the op- 
erators are represented as infinite-dimensional discrete 
matrices. We assume that the dependence S{t) is suffi- 
ciently smooth, so that it induces only short-range tran- 
sitions between quasi-energies [lo| . Then the operator in 
© rapidly approaches 1 away from the Fermi level, and 
the determinant converges for all values of A. Thus we 
conclude that x(u) is a single- valued analytic function of 
u = e*^ for any u e C\{0}. 

4. Determinant formula in z-representation. For our 
discussion, it will be convenient to rewrite ([3]) in a differ- 
ent form. We parametrize the 2AI x 2M unitary S'-matrix 
by two complex 2M x M matrices z and z as 



= iz\S) , 



(5) 



subject to the constraints z^z = z^z = 1m and z^z = 
Om- These constraints guarantee the unitarity of the S'- 
matrix. 



Substituting (O into ([3|) and using the explicit form of 
the charge operator P3, we immediately obtain 



x(A) — det (e" 



-iXnp ^tgi-^"* 



(6) 



Notice that the determinant size in the channel space is 
M X M in jH) versus 2M x 2M in ([3|. 

Furthermore, x(A) has a symmetry with respect 
to the right gauge transformation z h-^ zUuit) 
for any time-dependent unitary M x M matrix 
UM(t). Under this transformation, x(A) gets multi- 
plied by det (^C/Me-*^"^C/|^e*^"^) = exp[iXTT (np ~ 

UMnFUlf)] = exp iX J ^ tr {UMidtUlj) . Due to the 

periodicity of UM{t), the integral in the last expression 
is an integer number, and it affects only the overall shift 
of the total pumped charge. We conclude, that, up to this 
integer, the FCS depends only on z modulo U{M) gauge 
transformations. The latter space may be described as 
U{2M)/[U{M)xU{M)] and contains 2A'P real param- 
eters [instead of 4M^ real parameters in 5 G U{2M)]. 
It may be convenient to parameterize it with N{t) = 
2zz^ — 1, a 2M X 2M traceless Hermitian matrix obeying 
the constraint N'^{t) = I2M (this matrix was introduced 
in Ref. Physically, the gauge invariance described 

above corresponds to the independence of FCS on the 
scattering separately in the right outgoing states . 

In addition to gauge rotations, ([6]) is also insensitive 
to global (time-independent) left rotations z t—f U2MZ, 
which correspond to global rotations of N(t) [l3| . 

5. Derivation of the main result. To derive our main 
result we rewrite ([5]) as 

X{u) = det (1 + (m-1 - l)nF) {1 + (u - l)n'p) , (7) 

where Up = z^npz, and we have used n'j^ — up. 

The first factor in ([7]) does not involve any informa- 
tion about the time-dependent scattering. Its only role 
is to provide a reference point: the position of the Fermi 
level. The dependence on the evolution of S'-matrix en- 
ters through Up , which has a M x M matrix structure and 
tends to Om and to 1m at infinite positive and infinite 
negative energies, respectively. 

We can show that fT]) is fully determined by the spec- 
trum of rip, up to an overall constant charge transfer. 
Indeed, we can rewrite Q as 



X(A) = det 

= det 
„*ATr 



(e-*^"^e*^"^) (e-*^"'F[l + (u - l)n'p 
|^g-iA«FgiA«j.^ det (e-^-^"'F[l + iu- l)n^]) 
("J-"^)det (e-*^"^[l-|- (u- l)n|,]) , (8) 



where the determinant of the product can be written as 
the product of the determinants, since both operators 
behave properly (tend to 1m) at infinite energies, and 
the last transformation is justified, since rip — np rapidly 
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tends to zero at infinite energies. The trace in the first 
factor of the last expression equals the total transferred 
charge (generally non-integer) Q, 

Tt {n'p ~ np) ^ Q . (9) 

The second factor in ([8]) now depends only on the spec- 
trum of Up. 

For a periodic time dependence of S(t), the operator 
rip may be represented in the quasi-energy basis as a 
discrete infinite Hermitian matrix. It tends to at large 
positive energies and to 1 at large negative energies. Its 
spectrum is discrete and takes real values between and 
1, with possible accumulation points at and 1 [13]. If 
we denote the eigenvalues of Up as n^, Eq.([8|) implies the 
following expression for the generating function: 

X{u) = e^'Q n + • (10) 

a 

Notice that while this formula involves A = — ilnu, it 
is a single-valued function of u on C\{0}. As explained 
above, this expression is valid up to a constant integer 
charge transfer. 

It is obvious from ([7]) that the spectrum of rip is in 
one-to-one correspondence with the positions of zeroes 
of x(w), 

Uo, = l-{nirK (11) 

Since the spectrum of rip is real and lies between and 
1, we conclude that the zeroes of x('") must all lie on the 
negative real axis. This argument finishes the demon- 
stration of the main result of this paper. 

6. Example of forbidden charge transfers. We can il- 
lustrate our result with an example of a charge transfer 
which cannot be realized in our scattering system. Con- 
sider the following particular example of FCS: 

x{u) = l- 2F + F{u + u-^), (12) 

where < F < 1/2. The corresponding non- vanishing 
probabilities p±i = F , po — 1 ~ 2F are normalized and 
non-negative. However, the roots of (|12p are real only for 
F < 1/4. Therefore, the statistics ^ with 1/4 < F < 
1/2 cannot be produced by any evolution of scattering 
matrix at zero temperature. 

7. Allowed charge transfers. While the condition of 
real negative zeroes of x(u) is a necessary condition for an 
allowed FCS, it is also apparently a sufficient condition. 
Indeed, for any finite set of negative real Uq, the statistics 

Xiu) = n ^ (13) 

a— 1 

(which is identical to (|10|1 with Ni = Q — "a) '^^^^ 
be trivially realized in a system with N2 + 1 channels by 
using single-particle-transfer pulses as described in Ref. 4 



(one channel per each plus one channel for the overall 
integer transfer A'^i). 

Moreover, this FCS may also be arbitrarily closely ap- 
proximated in a single- channel system by sending indi- 
vidual well-separated pulses for each of the required ze- 
roes Ua. Indeed, in a single-channel system, the FCS is 
determined by the time evolution of the vector Tr Na 
on the two-dimensional sphere, as shown in Ref. 12). 
Motion of this vector along a circular trajectory of a 
given area (on the sphere) is equivalent to the "bias- 
voltage" problem with a fixed channel transparency de- 
termined by the area enclosed by the contour. This "bias- 
voltage" problem has been considered in Ref. where 
it has been shown that an arbitrarily sharp Lorentzian 
voltage pulse produces the elementary charge transfer 
x(u) — {u — Ua)/{l — Ua). By supcrimposing such pulscs, 
well separated in time, along circles of different areas, we 
can approximate the required statistics (jlSp arbitrarily 
close. 

Furthermore, since we can make this approximation for 
any finite number of roots N2 , and the possible infinite 
sets of roots have accumulation points only at Ua — > 
and Ua — » —00, we can also approximate arbitrarily close 
any statistics (fTO|) with an infinite set of roots, provided 
they converge sufficiently rapidly to — > and — * 1 
(for the convergence of (fTO|) it is sufficient to require that 
converge not slower than \a\~^ with > 1 at their 
accumulation points). We leave it as a mathematical 
problem to determine the precise requirements on the 
convergence and the conditions under which the generat- 
ing function (fTO| may be reproduced exactly by a suitable 
evolution of the S'-matrix. 

8. Bias voltage case. With this result in mind, we 
would like to comment on the existing results for the 
problem with a restricted evolution of the scattering ma- 
trix, namely the "bias- voltage" case \d\ . In this restricted 
problem, the transparencies of the channels are fixed, and 
the evolution of the 5-matrix is determined by the single 
parameter (/)(t) = J V{t)dt. The time dependence of z in 
this setup is given by z = e"^'-*'^-'zo, where zq is a fixed 
2M X M matrix. As follows from the results of Ref. 6, for 
this restricted evolution of the scattering matrix, there is 
an additional constraint on the positions of Ua. Namely, 
there are two t^pes of roots Uq : ordinary ( "typical" in no- 
tation of Ref. 6) and anomalous. Ordinary roots come in 
inversion-symmetric pairs (ua, 1/uq). In addition, there 
are M anomalous roots each having the same multiplicity 
\W\, where W is the winding number of 4'{t) (we assume 
that it is integer). They are located at either —gi/ (l — gi) 
or at —(1 — gi)/gi, depending on whether W is positive 
or negative. Here gi are the channel transparencies given 
by the eigenvalues of ZqPsZq. 

We can easily rederive this result within our formalism 
with a procedure analogous to that in Ref. [Q. For sim- 
plicity we consider only one channel with transparency 
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g. Then 



If = (1 - 9)nF + gn%, 



(14) 



where rip — e^'^'^'upe^''' . The symmetry of the spectrmii 
of n^p (at zero temperature) can be demonstrated with 
the use of the algebra of four operators, 



Q+ = np 
Qs = i[np,n'^] = 



iQ-Q-{ 



Q- = np 



.p, (15) 

1-Q+. 



The operators Q+i Q-i Q3 anticommute with each 
other, and the operator C commutes with them (this 
algebra rehes only on n^? and rip being projectors). Ob- 
viously C also commutes with 



1 



.9 0- 



1 



± 



-.9(1 -.9)C 
(16) 

and they can be diagonalized simultaneously [the last 
equality in pB]) should be understood as a relation be- 
tween eigenvalues, with different sign choices for different 
eigenvectors]. Let ^' be a common eigenvector of C and 
rip with the eigenvalues Ca and n^, respectively. The 
vector Qs^, if nonzero, has the eigenvalue 1 — n^, which 
produces [according to (fTTj) ] a pair of ordinary roots (ua , 
1/ua). If Qs'^ — 0, then one easily proves that either 
Q-^f or Q+\E' is zero. If Q-* = 0, then Ca = 0, and 
this corresponds to equal or 1, without any contri- 
bution to FCS. Finally, the zero modes Q+'^f — produce 
an anomalous contribution with Ca = 1. 

One can prove that, for a phase winding VF, there are 
exactly zero modes of Q+. Assuming a non- negative 
W (without loss of generality) , zero modes of Q+ must be 
simultaneously eigenstates of and of rip with eigen- 
values 1 and 0, respectively. Such eigenstates are eas- 
ily constructed explicitly by using the decomposition 



^2[Wut+4,+ {t)+4,^{t)]^ .^^j^gj.g 0^ g^j.g ^.j^g positive- 
and negative- frequency parts, and uj is the driving fre- 
quency. One finds that the space of zero modes is then 
spanned by Vfc = e-''"^*-">'-^*^ for fc = 1, 2, . . . , M^. The 
corresponding eigenvalue does not have a pair, but is 
IF-fold degenerate — 1 ~ g. Similarly, for W < 0, one 



and n 







finds \W\ zero modes of with np 
which produces = g. 

In terms of the generating function, the separation into 
the inversion-symmetric (ordinary) and anomalous (de- 
pending only on the average charge transfer) parts may 
be written as (for > 0) 

X{u) = [1 + giu - 1)]^ Xinv(u), Xinv{u^^) = Xinv(u) 



With some algebraic manipulations, one can also derive 
an explicitly inversion-symmetric formula for Xu 



Xinv(M)=det 1 + gil - g){u + ■ 



2)(l-nF)4 



(this expression also assumes W > 0). This formula may 
be viewed as a compact form of the product over eigen- 
values obtained in Ref. @. Note that this separation of the 
two contributions is specific to the "bias-voltage" case at 
zero temperature. 

9. Conclusion. We have derived constraints on the 
charge-transfer statistics between two conductors in the 
limit of instant scattering and at zero temperature. Our 
findings generalize the results obtained previously for the 
case of a fixed scatterer with a time-dependent voltage 
The allowed statistics is characterized by negative 
real zeroes of the characteristic function x{u)- While 
we have not performed a similar analysis of the problem 
at finite temperature, we conjecture that the singulari- 
ties of logx(u) remain restricted to the negative real axis 
even at finite temperature, but develop a cut instead of 
a discrete set of branching points. The available sim- 
plest examples of the charge-transfer statistics at finite 
temperature support this conjecture. 

Similarly to Ref. i6|, our result pO)) can be interpreted 
as a decomposition of the charge transfer into indepen- 
dent tunneling events - binomial processes. The effective 
channel transparencies of these events are the eigenvalues 



of; 



depending on the time evolution of the scattering 



matrix. This interpretation suggests that the physical 
origin of this result lies in the assumed absence of inter- 
actions between electrons and that it may break down 
once interaction is taken into account 15 1. 
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